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ABSTRACT 

Kinetic plasma theory is used to generate synthetic spacecraft data to analyze and interpret 
the compressible fluctuations in the inertial range of solar wind turbulence. The kinetic coun- 
terparts of the three familiar linear MHD wave modes — the fast, Alfven, and slow waves — are 
Q-t identified and the properties of the density-parallel magnetic field correlation for these kinetic 

wave modes is presented. The construction of synthetic spacecraft data, based on the quasi-linear 
premise — that some characteristics of magnetized plasma turbulence can be usefully modeled as 
a collection of randomly phased, linear wave modes — is described in detail. Theoretical predic- 
tions of the density-parallel magnetic field correlation based on MHD and Vlasov-Maxwell linear 
cigenfunctions are presented and compared to the observational determination of this correlation 
based on 10 years of Wind spacecraft data. It is demonstrated that MHD theory is inadequate to 
describe the compressible turbulent fluctuations and that the observed density-parallel magnetic 
field correlation is consistent with a statistically negligible kinetic fast wave energy contribution 
for the large sample used in this study. A model of the solar wind inertial range fluctuations is 
proposed comprised of a mixture of a critically balanced distribution of incompressible Alfvcnic 
fluctuations and a critically balanced or more anisotropic than critical balance distribution of 
compressible slow wave fluctuations. These results imply that there is little or no transfer of 
large scale turbulent energy through the inertial range down to whistler waves at small scales. 

Subject headings: turbulence — solar wind 



1. Introduction 

Despite more th an forty years of direct spacecraft measurements of turbulence in the near-Earth solar 



wind (|Colemanlll968l ). our understanding of turbulence in a magnetized plasma remains incomplete. One of 
the primary goals is to understand the role of the turbulence in mediating the transfer of energy from large 
to small scales. Within the turbulent inertial range of scales, corresponding to spacecraft-frame frequencies 
of 10~ 4 Hz < f sc < 1 Hz or length scales 10 6 km > A > 10 2 km, the fluctuations involved in this energy 
transfer are a mixture of co mpressible and incompressible fluctuations, w ith around 90% of the energy in the 



incompressible component ( Tu fc Marsch 1995 ; Bruno fc Carbonej 2005 ) . These incompressible fluctuations 
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have been identifed as Alfven waves ( Belcher fe Davisl 1971 ). but the nature of the compressible component 
remains uncertain. 

The compressible turbulent fluctuations have often been interpreted as a combination of magneto acoustic 
(fast MHD) waves and pressure-balanced structures (|Tu fe Marschlll995t iBruno fe Carbond 120051 ) . Early 
studies of thermal and magnetic pressure fluctuations in the solar wind found an anti-correlation of the 
thermal pressure and magnetic pressure at timesc ales of 1 h, corresponding to an interval of constant total 
pressure, or a pressure- balanced structu re (PBS) (Burlagalll968 ; Burlaga fe Qgilvie 1970 ). Further studies 
found evidence of PBSs out to 24 AU (jBurlaga et alJll990| ). Related investigations discovered a similar 
anti-correlation between the density n and magnetic field magnitude B from 0.3 AU to 18 A U on timescale s 
ranging from several hours to 1.8 minutes ( Vellante fc Lazarus! 1987 : Roberts et al. 1987al fbl: Roberts 1990h . 
Theoretical studies of compressive MHD fluctuations in the low-Mach number, high-/? limit interpreted 
these anti-co rrelated density-magnetic field strength observati ons as nonpropagating "pseudosound" density 
fluctuations ([Montgomery et al.lll987t iMatthaeus et al.lll99ll ). Later, a more comprehensive observational 
investigation confirmed the general density-magnetic field magnitude anti-correlation, but also identified a 
few positively correlated intervals consistent with the magnetosonic (fast MHD) wave (|Tu fc Marschlll994 h 
Analy sis of Ulysses observations found evidence for PBSs at inertial ran ge scales in the high latitude solar 
wind ( McComas et al.lB-995 ; Reisenfeld et al.lll999 ; Bavassano et alj|2004 ). Studies of the electron density up 
to / = 2.5 Hz also found pressure balanced structures but interpreted these as ion acoustic (slow MHD) waves 
and r ecognized that PBSs are simply the ion acoustic (slo w MHD) wave i n the perpendicular wavevector 
limit (jKellogg fc Horburvll2005l ). a fact previously noted bv iTu fc Marschl (|1994h . Recently, measurements 
of the anti-correlation between electron density an d magnetic field strength indicated the existence of PBSs 
over timescales ranging from 1000 s down to 10 s (|Yao et al.ll201lh . 



An important consideration in the study of the compressible fluctuations of solar wind turbulence is 
the fact that the mean free path in the solar wind plasma is about 1 AU, so the dynamics over the entire 
inertial range is weakly collisional. The implications of this fact have not been seriously addressed in any 
of the aforementioned studies of compressible fluctuations in solar wind turbulence. The MHD description 
is rigorously valid only in the limit of strong collisionality, so a kinetic description is formally required to 
describe t he inertial range turbulence. In the limit fey <C predicted by anisotropic MHD turbulence 
thcoric fEI (jColdreich fc Sridharl[l99i iBoldvrevlEooil . it has been demonstrated that even in the weakly 
collisional limit, the turbulent dynamics of the Alfven waves decouples from the comp ressible fluctuations 
and is rigorously described by the equations of reduced MHD (jSchekochihin et al.ll2009l ). The compressible 
fast and slow wave modes, on the other hand, require a kinetic description to resolve both the wave dynamics 
and the collisionless kinetic damping mechanisms. The study presented here is the first to examine the 
properties of the compressible fluctuations in the turbulent solar wind using Vlasov-Maxwcll kinetic theory. 
Specifically, we use the predicted correlation between the density fluctuations and parallel magnetic field 
fluctuations to determine the nature of the compressible fluctuations in the solar wind. 

In Sj2l we explore the connection between the familiar linear wave modes in MHD and the corresponding 
kinetic wave modes in Vlasov-Maxwcll kinetic theory and demonstrate that weakly collisional conditions do 
not change the qualitative properties of the density-parallel magnetic field correlation. In we discuss the 
quasi-linear premise upon which the method of synthetic spacecraft data is based and describe in detail the 
procedure for generating synthetic spacecraft data. The synthetic spacecraft data predictions of the density- 
parallel magnetic field correlation based on linear MHD and Vlasov-Maxwell eigenfunctions is presented in 



1 Parallel and perpendicular are denned with respect to the direction of the local mean magnetic field. 
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Sj4j A comparison of the synthetic spacecraft data predictions to the observational determination of the 
density-parallel magnetic field correlation is presented in $5J showing a statistically negligible fast wave 
energy contribution to the compressible fluctuations. The implications of this finding are discussed before 
summarizing the findings of this investigation in SJ5] 



2. Collisional vs. Collisionless Dynamics of Compressible Fluctuations 



The inertial range of solar wind turbulence is observed to be a mixture of inco mpressible and com- 
pressible motions, with a round 90% of the energy due to the incompressible component (|Tu fe Marscblll995 ; 
Bruno fc Carbone] 2005 ). If these fluctuations are interpreted as some mixture of the three MHD linear 
wave modes, then Alfven waves are responsible for the incompressible component, while slow and fast MHD 
waves make up the compressible component. In the limit of large scales compared to the thermal ion Larmor 
radius, kpi <C 1, these modes may be distinguished by the correlation between the density and parallel 
magnetic field fluctuations: fast waves are positively correlated, slow waves are negatively correlated, and 
the density and parallel magnetic field fluctuations are both zero for Alfven waves. As the wave amplitude 
is increased to nonlinear levels, even in the limit that they form discontinuities or shocks, these qualita- 
tive properties persist, correspon ding to tangential and rotational discontinuities or fast and slow shocks 
(jBaumiohann fc Treumannlll996T ). In this section, we will explore the properties of the kinetic counterparts 
to the fast and slow MHD wave modes in the inertial range using Vlasov-Maxwell kinetic theory. 



2.1. The Inertial Range Limit of the Compressible Linear Wave Modes 



As derived in Appendix [XI the normalized compressible MHD linear dispersion relation depends on 
only two parameters, u)/{kvA) = ujmhd{P,9): the plasma beta, j3 = c 2 s /v\, where the sound speed is 
c s = v 7Po/po and the Alfven velocity is va = Bq / y/Airpo; and the angle 9 between the local mean magnetic 
field Bo and the direction of the wavevector. 

To establish precisely the connection between the three linear MHD wave modes and their kinetic 
counterparts in the linear Vlasov-Maxwell system, we specify a fully-ionized proton and electron plasma 
with isotropic Maxwellian velocity distributions and a realistic mass ratio rrii/m e = 1836. In general, the 
linear Vlasov-Maxwell dispersion relation depends on five parameters: the ion plasma beta = 87rnjTj/ Bq, 
which is the ratio of the ion thermal pressure to the magnetic pressur^f], the normalized wavenumber kpi, 
the angle 9 between the local mean magnetic field Bo and the direction of the wavevector, the ion to electron 
temperature ratio Ti/T e , and the ratio of ion thermal velo city to the speed of light yule. The solutio n 
may then be expressed as u = uivAfiPi, kpi,9,Ti/T e ,vu/c) (Stix Il992t lQuataertlfl995 iHowes et allbood ). 
To connect to the single fluid theory of MHD, we take equal ion and electron temperatures, Ti/T e = 1. 
In the inertial range limit, kpi <C 1, and for the non-relativistic conditions appropriate to the solar wind, 
vu/c <C 1, the normalized linear Vlasov-Maxwell dispersion relation simplifies to Lu/(kvA) = wyA/ (A, 
Since (3 = (3i + (3 e = (3i(l + T e /Ti), direct quantitative comparison betwee n the solutio ns of the MHD and 
Vlasov-Maxwell linear dispersion relations is possible by choosing (i = 2/3, I Howes! 2009 ). 



The Vlasov-Maxwell solutions presented in this section use the parameters kpi = 2.2 x 10 4 , — = 10 4 , 
and Ti/T e = 1, and employ a realistic mass ratio mi/m e = 1836. The complex eigenfrequencies are solved 



2 In this study, the Boltzmann constant is absorbed into the temperature to yield temperature in units of energy. 
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numerically ( Quataert 1998 ; Howes et al. 20061) with Bessel function sums evaluated to 100 terms to ensure 
accurate results. 



2.2. Connection between MHD and Vlasov- Maxwell Linear Wave Modes 



Previous studies have co mpared the properties of lin ear Vlasov-Max well wave mo des with the linear 
modes from two-fluid theory (jKrauss-Varban et al.l Il994f) and Hall MHD (jHowesI 120091 ) . Here we restrict 
ourselves to establishing the connection between MHD wave modes and Vlasov-Maxwell wave modes in 
the inertial range limit, kpi <C 1. Due to the reduced parameter space of the Vlasov-Maxwell system in the 
incrtial range limit, the comparison between the three linear MHD wave modes and their kinetic counterparts 
is concisely expressed by the use of normalized Freidrichs diagrams, polar plots of the wave phase velocity 
normalized to the Alfven speed ui/{kvA) vs. polar angle 9, as shown in Figure [T] for four values of ft. 



As discussed by iKrauss-Varban et alj (| 1994T) , the identification of corresponding wave modes is com- 
plicated by the existence of a branch cut in the complex solution space of the Vlasov-Maxwell dispersion 
relation. At inertial range scales kpt -C 1, this branch cut exists only at small angles 9; for wavevectors in 
the perpendicular limit, the kinetic fast, Alfven, and slow modes are always easily distinguished. Therefore, 
we adopt the strategy that we label the modes according to their properties in the perpendicular wavevector 
limit, and retain those labels for the linear dispersion relation solutions as the angle is decreased from 9 = 90° 
to 9 = 0°. Further discussion of the properties of the kinetic fast and slow wave modes at small angles 9 is 
presented in Appendix [Bl 

Results of the MHD to Vlasov-Maxwell comparison are presented in Figure [T] In the left column are 
the normalized Freidrichs diagrams for the MHD results for the fast wave (red), Alfven wave (black), and 
slow wave (blue). In the center column are the Vlasov-Maxwell results for the kinetic fast wave (red), Alfven 
wave (black) , and kinetic slow wave (blue) . Note that the strongly damped Vlasov-Maxwell modes (defined 
by —"f/uj > 0.25, where uj + i"f is the complex Vlasov-Maxwell cigcnfrequcncy) arc given by dashed lines. In 
the right column are the paths of the Vlasov-Maxwell solutions in complex uj space, running from the nearly 
parallel limit (crosses) to the nearly perpendicular limit (open circles). The region below the grey dashed 
line is — 7/w > 0.25, indicating strong collisionlcss damping. Rows present different values for ion plasma 
beta, ft = 0.1,1.25,1.58,10.0. 

For most values of ft and 9 presented in the leftmost two columns of Figure [IJ the correspondence 
between the MHD and kinetic wave modes is clear. It is worth noting, however, several distinctions between 
the fluid and kinetic behavior. First, the kinetic slow wave has a greater phase speed than the Alfven wave 
for ft > 1. In fact, for sufficiently parallel wavevectors and ft ~ 1, the kinetic slow mode also has a greater 
phase velocity than the fast mode. Second, the typical magnctic-to-acoustic mode conversion that occurs 
for the MHD fast and slow modes at /3 = 1 is replaced by a conversion between the compressible kinetic 
roots at ft ~ 1.3. This conversion leads to damping of the fast waves at small angles 9 in the large ft limit. 
This general identification of the kinetic fast and slow wave modes is used in the construction of synthetic 
spacecraft data, as described in $3l 
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Fig. 1. — Fredrichs Diagrams of normalized phase speed uj/(kvA) for MHD (left column) and Vlasov- 
Maxwell (middle column) modes for kpi = 2.2 x 10 -4 and values of f3 = 2/3^ = 0.2, 2.5, 3.16, and 20. Parallel 
wavevectors, 9 = 0° with respect to the local mean magnetic field, lie along the the vertical axis, while 
perpendicular wavevectors, 8 = 90°, lie along the horizontal axis. Dashing indicates — 7/w > 0.25, which 
implies heavy collisionless damping of that wave mode. On the right are the paths of the Vlasov-Maxwell 
solutions in complex to space, running from the nearly parallel limit (crosses) to the nearly perpendicular 
limit (open circles), with j/uj = —0.25 plotted as a dashed grey line. 
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Fig. 2. — The normalized correlation by Fourier mode, Ck(5n, SBu), for the kinetic fast (left), Alfven 
(center), and slow (right) modes over the plane (9, Pi). A positive correlation is found for the kinetic fast 
modes, and a generally negative correlation for the slow modes. 

2.3. Density-Parallel Magnetic Field Correlation for Kinetic Fast and Slow Modes 

With a clearly defined identification of the kinetic fast and slow modes complete, we may now calculate 
the density-parallel magnetic field correlation for the kinetic fast and slow modes. This will enable us to 
verify whether the general qualitative properties of this correlation — that fast waves are positively correlated, 
and slow waves negatively correlated — remain unchanged for collisionless conditions. The density fluctuation 
8n and parallel magnetic field fluctuation 5B\\ may be calculated from the complex solution of the linear 
Vlasov-Maxwell dispersion relation. For a chosen wavevector k, we obtain complex Fourier coefficients 5n(k) 
and 8Bu (k), and we define the normalized correlation by Fourier mode, Ck, given by 

<w^-(i£»)- 

Note here that the Fourier coefficients of both density and parallel magnetic field fluctuations satisfy the 
reality condition, e.g., <5n(k) = Sn*(— k). 

In Figure [3 Ck(8n, 5B\\) is plotted for the kinetic fast (left), Alfven (center), and kinetic slow (right) 
modes as a function of the parameters (9, Pi). The remaining parameters for the Vlasov-Maxwell eigen- 
function solutions are kpi = 2 x 10 -4 , Ti/T e = 1, and Vtjc = 10~ 4 . The kinetic fast mode (left) always 
has a positive correlation, with Ck(Sn, 5B\\) > 0.9 over most of the (9, Pi) plane. The kinetic slow mode 
(right) has a negative correlation with Cu(5n, SB^ ) < —0.9 over most of the {9, Pi) plane; the correlation 
becomes slightly positive for 9 < 30° and Pi > 1.3. The Alfven mode (center), presented for completeness, 
has a more complicated behavior, but it is worthwhile noting that the amplitudes of <5n(k) and <5S||(k) are 
both very small (compared to the characteristic amplitudes for either of the compressible waves) at incrtial 
range scales, so this correlation for Alfven waves is likely to be unmeasurable in the solar wind. In general, 
these results confirm that the general qualitative properties of the density-parallel magnetic field correlation 
remain unchanged for the kinetic fast, Alfven, and slow modes in a weakly collisional plasma. 



3. Constructing Synthetic Spacecraft Data 

The primary methodology employed in this study of the compressible fluctuations in the solar wind 
inertial range is the construction of synthetic spacecraft data for direct comparison to actual single-point 
spacecraft measurements. For the present study of the correlation between the density and parallel magnetic 
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ficld fluctuations, this technique enables us to determine the characteristics of the correlation in the presence 
of a turbulent spectrum of wave modes. This method of analysis is based on the quasi-linear premise that 
some characteristics of magnetized plasma turbulence can be usefully modeled as a collection of randomly 
phased, linear wave modes; the justification for this premise is discussed in ij3.ll The result from Q2.1\ — 
that, in the inertial range, the dynamics of both MHD and Vlasov-Maxwell plasmas depend only on two 
parameters, the ion plasma beta /3j and the wavevector angle 9 — significantly simplifies the construction 
of synthetic plasma data. Upon adoption of the quasi-linear premise and specification of the plasma /3j, 
one need only choose the partitioning of power among the contributing linear wave modes (fast, Alfven, or 
slow) and the wavevector distribution of power (isotropic or critically balanced) for each of those modes, as 
described in ^3.21 



3.1. The Quasi-linear Premise: Modeling Turbulence as a Spectrum of Linear Wave Modes 



The theoretical investigation of turbulence usingsynthetic data is based on a concept that we denote 
the quasi-linear premise. The quasi-linear premise □ states that some properties of magnetized plasma 
turbulence can be understood by modeling the turbulence as a collection of randomly phased, linear waves. 
In this picture, the nonlinear turbulent interactions serve to transfer energy from one linear wave mode 
to another — thus, the picture is quasi-linear. The mathematical properties of the equations that describe 
turbulence in a magnetized plasma, in conjunction with a phenomonological understanding of the properties 
of the turbulence, provide the motivation for this quasi-linear approach. This concept can be most easily 
explained using the following example of turbulence in an incompressible MHD plasma. 



The Elsasser form of the ideal incompressible MHD equations (|Elsasserlll950l ) is given by 



^ T (v A .V)z ± = -(z^.V)z ± -^, 
dt p a 



(2) 



where the magnetic field has been decomposed into its equilibrium and fluctuating parts B = Bo + <5B, 
the Alfven velocity due to the equilibrium magnetic field Bo is given by va = 'Bo/y/iirpo, and z^ = 
Sv ± SH/y/Airpo are the Elsasser fields describing the velocity and magnetic field behavior of waves traveling 
down (up) the mean magnetic field. The second term on the left-hand side of equation (0) represents the 
linear propagation of the Elsasser fields along the mean magnetic field at the Alfven speed, while the terms on 
the right-hand side represent the nonlinear interactions between upward and downward propagating waves, 
where the pressure gradient term ensures incomprcssibility of the fluctuations. 



The theory of strong incompressible MHD turbulence (jGoldreich fc Sridharlll995[ lBoldvrevll2006l ) sug- 
gests that the turbulent fluctuations at small scales become anisotropic, where the nonlinear cascade of energy 
generates turbulent fluctuations with smaller scales in the perpendicular direction than in the parallel di- 
rection, fc|| <C k±. This inherent an isotropy of magnetized plasma turbulence has long been recognized from 



early studies in laboratory plasmas ( Robinson fc Rusbridg 



1971: Zweben et al.lll979: 



Montgomery fc Turner 



19811 ) and in early numerical simulations (jShebalin et al.l ll983l). It has been conjectured that strong turbu- 
lence in incompressible MHD plasmas maintains a state of critical balan ce between the linear timescale for 
Alfve n waves and the nonlinear timescale of turbulent energy transfer (|Higdonlll984t iGoldreich &: Sridhar 



19951 ). There exists significant evidence consistent with the predictions of critical balance from analysis 



3 Note that the quasi-linear premise is not the same as, and does not require, the quasi-linear approximation, a rigorous 
mathematical procedure that requires weak nonlinear interactions such that perturbation theory can be applied. 
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of numerical simulations (Cho & Vishniac 2000; 


Maron & Goldreich 2001; Cho & Lazarian 2003) and solar 


wind observations 


Horburv et al. 20081 Podesta 


2009: Wicks et al. 201(A Luo & Wu 2010: Chen et al. 


2011: 


Forman et al. 


2011 


). In a state of strong turbulence, critical balance implies that the linear term (va ■ V)z ± 



and nonlinear term (z^ • V)z ± in equation ^ arc of the same ordeiQ It is this property that motivates the 
adoption of the quasi-linear premise, as discussed below. 

In the absence of the nonlinear terms (setting the first term on the right-hand side of equation to zero), 
the behavior of the plasma is entirely determined by the linear term. If the right-hand side of the equation is 



of the plasma to the imposed perturbation 



Montgo mery fc Matthaeus 



1997 



Galtier et alJl200 



E3 



1995 



Lithwick fc Goldreich 



of weak turbulence i 


Sridhar & Goldreich 


1994: 


Goldreich & Sridhar 


1997: Ntr & Bhattacharjee 



200 



J| the nonlinear terms on the right-hand of equation ([2]) 



are indeed a small perturbation to the linear system, representing the nonlinear transfer of energy between 
the linear wave modes. Perturbation theory may be applied to the study of the turbulent dynamics in this 
limit, so the quasi-linear premise is clearly valid for the case of weak turbulence. 

For strong turbulence, the condition of critical balance implies that the energy in a particular linear 
wave mode may be transferred nonlinearly to other modes on the timescale of its linear wave period. But 
since the linear and nonlinear terms are of the same order in critical balance, the linear term still contributes 
significantly to the instantaneous response of the plasma, even in the presence of strong nonlinearity. There- 
fore, the fluctuations in a strongly turbulent magnetized plasma are expected to retain at least some of the 
properties of the linear wave modes. In particular, for a turbulent fluctuation with a given wavevector, the 
amplitude and phase relationships between different components of that fluctuation are likely to be related 
to linear cigenfunctions of the characteristic plasma wave modes. In the construction of synthetic plasma 
turbulence data, a spectrum of randomly phased linear wave modes can be specified, with the amplitude of 
each of the linear modes adjusted to satisfy a chosen observational constraint, such as the turbulent magnetic 
energy spectrum (the second-order moment). By adopting the quasi- linear premise, the properties of the 
synthetic turbulence data may then be compared directly to spacecraft measurements to explore the nature 
of turbulent fluctuations. 

The third- and higher-order moments of the turbulence, on the other hand, are clearly not described by 
this simplified quasi-linear approach. Such higher order statistics depend critically on the phase relationships 
between different linear wave modes, and these phase relationships are determined by the nonlinear interac- 
tions responsible for the turbulent cascade of energy from large to small scales. For a collection of randomly 
phased linear waves, such higher order statistics of synthetic data constructed using the quasi-linear premise 
will average to zero, yielding no useful information. 

Although our illustration of the application of the quasi-linear premise above specifics the case of turbu- 
lence in an incompressible MHD plasma, the necessary general properties of the linear and nonlinear terms, 
as well as the inherent anisotropy of magnet ized plasma turbu l ence, continue to hold fo r less restricted 
plasma conditions, including kinetic plasmas ( Howes et al. 2006 : Schekochihin et al. 20091 ). In particular, 
the arguments for the importance of the linear physics even in a strongly turbulent plasma also hold for 
the linear collisionless damping rates, providing the foundation for simple models of the turbulent energy 
cascade, encompassing both the inertial and dissipation ranges (jHowes et al.ll2008ai 120111 ). 



4 Note that, in the case of incompressible hydrodynamic turbulence (Euler or Navicr-Stokes), the absence of a linear term 
prohibits the possibility of a quasi-linear approach. 
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Besides the feasibility arguments for the validity of the quasi-linear premise outlined above, we give here 
no a priori proof for its validity in strongly nonlinear plasma turbulence. Nonlinear simulations of plasma 
turbulence and observational studies of solar wind turbulence provide two avenues for testing the validity of 
the premise — at present, there exist arguments in the literature both for a nd against its validity. Analysis of 



nonli near numerical simulations of plasma turbule nce using gyrokinetic s (jHowes et al.ll2008bt iHowes et al 



iity 

m 



2011 ) and both Hall MHD and Landau fluid theory Jrlunana et aL 2011 ). as well as observational analysis of 
multi-spacecraft data in the solar wind ( Sahraoui et alJl2010h . support the valid ity of the quasi-linear 



premise. 



In addition, given that the idea of critical balance in strong MHD turbulence (jGoldreich fe Sridharlll995l ) is 
essentially a quasi- linear concept — that the timescale of the non linear energy transfer remains of order the 



linear wave frequency— evidence in support of critical balance (ICho fe Vishniac 



200ll: ICho fe Lazarianlbooi iHorburv et ZII2OO8I: IPodestalliool IWicks et al 



2010: 



200d iMaron fe Goldreich 



Forman et al]|201lf) also 



i ndirectly supports the quasi -linear premise. In contrast, studies of 3D incomp ressible MHD simulations 
( Dmitruk fe Matthaeusl 20091) and 2D hybrid simulations ( Parashar et al. 2010|) of plasma tu rbulence, as 
well as an observational analysis of multi-spacecraft data in the solar wind (jNarita et al.ll201ll) . have called 
into question the validity of the quasi-lin ear premise. A revi ew of this supporting and conflicting evidence for 
the quasi- linear premise is presented in IHowes et al.l ( 2011 ). focusing in particular on questionable aspects 
of the conflicting studies that cast doubt on the validity of their conclusion that the quasi-linear premise is 
inapplicable to the case of strong plasma turbulence. 

It is also important to note that the utility of the quasi-linear premise for the study of plasma turbulence, 
however, may also be judged a posteriori by the insights gained from such an approach. 



3.2. Procedure for Constructing Synthetic Spacecraft Data 

Upon adopting the quasi-linear premise, the construction and analysis of synthetic spacecraft data 
requires three steps: 

1. Populate a synthetic plasma volume with a spectrum of linear wave modes with a chosen distribution 
of power in wavevector space. 

2. Sample the synthetic plasma volume at the position of a probe moving with respect to the plasma to 
generate reduced time series comparable to single-point spacecraft measurements. 

3. Perform the requisite analysis on the synthetic time series to compare to spacecraft data analysis, for 
example, the zero-lag cross-correlation of the density and parallel magnetic field fluctuations. 

For the study of the density and parallel magnetic field correlation of the solar wind compressible turbulent 
fluctuations, each of these steps is detailed below. 



3.2.1. Creating the Synthetic Turbulent Plasma 

As discussed in N2. 11 for the inertial range of solar wind turbulence, the appropriately normalized linear 
physics of both the MHD (collisional fluid) and Vlasov-Maxwell (collisionless kinetic) systems depends on a 
reduced set of two parameters: the ion plasma beta Pi and the wavevector angle 9. Therefore, a completely 
general synthetic turbulent plasma requires the specification of just three properties: (1) the plasma beta, (2) 
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the fraction of power in each of the possible linear wave modes, and (3) the distribution of power in wavevector 
space for each of these wave modes. Both observational constraints and phenomcnological models of plasma 
turbulence guide our choices for these properties. 

First, we discretize the three-dimensional wavevector space on a uniform grid of 32 3 points, where each 
wavevector component spans kipi £ [—4.8 x 10~ 2 ,4.8 x 10~ 2 ] with a minimum grid spacing fcoPi = 3x 10 -3 . 
Taking an equilibrium magnetic field B = B a z and specifying the ion plasma beta we may solve for 
the normalized linear frequencies uJ(k) and the linear eigenfunctions £B(k), <JE(k), <5U(k), <5n(k) for each 
wavevector k using our chosen plasma description (MHD or Vlasov-Maxwell) . Note that we adopt the 
convention uj > 0, so that the direction of the wave group velocity (in the case of Alfven waves, up or down 
the mean magnetic field) is determined by the wavevector. 

Next, we specify the fraction of the turbulent magnetic power for each Fourier component due to the 
combination of the fast, Alfven, and slow waves. After specifying the parititioning of turbulent power among 
the linear wave modes, we adjust the amplitudes of the Fourier coefficients of the linear wave modes so 
that the fluctuating magnetic power |<5B(fc)| 2 (due to the linear superposition of all of the contributing wave 
modes at a given wavenumber k) is consistent with the inertial range observational constraint that the one- 
dimensional magnetic energy spectrum scales as Eb <x A;" 5 / 3 . The amplitudes of the Fourier coefficients for 
the remaining fields <5E(k), <5U(k), and <5n(k) are specified by the eigenfunction solution of each of the linear 
wave modes. Random phases are also applied to each wave mode. 

Finally, the distribution of energy in wavevector space for each of the constituent wave modes must be 
specified to model the inherently anisotropic nature of magnetized plasma turbulence. Imbalance between 
the turbulent energy fluxes propagating up and down the magnetic field as well as anisotropy in the angular 
distribution of turbulent energy with respect to the mean field direction may be incorporated in this final 
step. For this initial synthetic data investigation of the compressible fluctuations in the inertial range, we 
always set the upward and downward propagating wave energy fluxes in balance. Numerical simulations of 
compressible MHD plasma turbulence suggest that fast wave energy is distributed isotropically while both 
the Alfven and slow wave energies obey a critically balanced distribution with energy con centrated mainly 



in m odes with ku < kV 3 k^/ 3 , where ko is the isotropic driving scale of the turbulence (|Cho fc Lazarian 



2003). Therefore, we allow for two possible wavevector distributions of energy: (a) an isotropic distribution, 



such that the distribution of power is independent of 9; and, (b) a simplified critically balanced distribution, 

1/3 2/3 

where all modes with fc|| > k k^_ are set to zero, where k^ is the minimum wavenumber of the simulation 
domain. 

Note that the adjustment of the fraction of turbulent power due to each wave mode is performed for 
each Fourier component, so that changes in the wavevector distribution of energy — e.g., from isotropic to 

1/3 2/3 

critically balanced by zeroing out all modes with kn > k fcj_ — do not affect the fluctuation amplitudes of 
the non-zero modes. 



3.2.2. Generating Synthetic Reduced Time Series 

Once the synthetic turbulent plasma has been completely specified in the Fourier domain, the turbulent 
fields may be computed at any position r and time t by summing over all contributing Fourier modes of all 
constituent linear waves, for example, 

m k 
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Here, the index m indicates the contributions from fast, Alfvcn, and slow modes, and each contributing 
Fourier mode k for each wave mode m is given a constant, random phase (f> m j, where j = 1, . . . ,n x n y n z . 
Note that the linear frequency uj m (k) for each constituent wave is a function of the wavevector k as well as 
the wave type m. 

Unfortunately, single-point satellite measurements do not provide full spatial information about the 
turbulence for comparison to the synthetic data. To mimic spacecraft measurements made as the supcr- 
Alfvenic solar wind streams past the satellite with velocity v, we sample the synthetic plasma volume at 
the position of a probe moving with velocity — v through the volume, r(t) = Tq — vt, where for simplicity 
we set Tq = 0. Sampling along the probe trajectory at an interval At generates single-point time series at 
times t n = nAt for each of the turbulent fluctuating fields, where n = 1, . . . N and the total time interval is 
therefore T = NAt. This reduced set of data is directly comparable to single-point spacecraft measurements. 
The time series of the parallel magnetic field fluctuation 8B\\ n = 8B\\(t n ), for example, is given by 

SB\\ n = 5Bu(r,f)| r= _ vtn =J2Yl SB \\™M e ~ i[{k ' V+Um)tn ~ M - W 

m k 



Note that the frequency u>' of the signal measured by the moving probe is Doppler shifted by the probe 
velocity, u>' = lu + k • v. Normalizing this Dopplcr-shiftcd frequency by kvA to obtain ul' = lo' /(Jcva), 
we find uJ' = uj + k • (v/u^), where k is the unit vector in the direction of the wave vector. The Taylor 
hypothesis — that the temporal fluctuations measured in the super- Alfvenic solar wind flow are dominated by 
spatial fluctua tions swept past the probe at the solar wind velocity is equivalent to the limit k- (v/va) ^> uJ 



iTavlor 1938 ). Because the ratio v/va — 10 in the super- Alfvenic solar wind flow (|Tu fc Marsch 



1995 



Bruno fc Carbong 2005 ). the Taylor hypothesis is frequently a useful simplification for studies of the non- 
dispersive linear wave modes of the inertial range. For the present study of the compressible fluctuations in 
the inertial range, we adopt the Taylor hypothesis that to' ~ k • v, so 

SB \\n = E E SB^e-^**-*^. (5) 

m k 

In £14.51 we test the effect that violation of the Taylor hypothesis has on the density-parallel magnetic field 
correlation C(Sn, SBh). 

This procedure of sampling the synthetic plasma volume is depicted schematically in Figure [3] Here 
the blue, green, and purple surface plots represent the spatial variation for a few of the contributing Fourier 
modes to the turbulent fields. The synthetic plasma volume is sampled at uniform time intervals along 
the probe trajectory (red line), generating single-point time series of the parallel magnetic field fluctuation 
SBu n (upper right) and the density fluctuation 5n n (lower right). These synthetic time series may then be 
analyzed using the same procedures as the actual spacecraft measurements. 



3.2.3. Analysis of Synthetic Spacecraft Data 

As discussed in ^2.3[ the correlation between the density fluctuation Sn and the parallel magnetic field 
fluctuation <5-B|| distinguishes fast from slow compressible modes in either the fluid MHD or kinetic Vlasov- 
Maxwell systems. The zero-lag cross-correlation of the two time series 8n n and 5B\\ n is given by 

C(Sn, £B|| ) = , =, (6) 



E n (5n n -5n) J^SB^- SB 
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Time 



Fig. 3. — Probe trajectory (red) through a synthetic turbulent plasma volume generated using a spectrum 
of linear eigenfunctions (blue, green, and purple surface plots on the left). Measurements of the plasma 
fluctuations along the probe trajectory yield synthetic single-point time series of the parallel magnetic field 
fluctuation 8Bu (upper right) and the density fluctuation 5n (lower right). 




where 5n and <5i3|| are the averages over the time interval T. Comparison of the theoretical predictions 
of C(5n, 5B\\) to the correlation from spacecraft measurements, assuming the validity of the quasi-linear 
premise, provides a means of constraining the nature of the compressible fluctuations in the solar wind. 



4. Theoretical Predictions of C(Sn, SB\\) Using Synthetic Spacecraft Data 

Our theoretical investigation of the compressible fluctuations in the inertial range of solar wind tur- 
bulence uses the method outlined in $3] to generate synthetic spacecraft data to provide predictions of the 
density-parallel magnetic field correlation C(<5n, 8B^ ) as a function of the following turbulent plasma proper- 
ties: the plasma /3,, the fraction of fast-to-total compressible wave power J 7 , and the angular distribution of 
the turbulent power in wavevector space. These theoretical predictions may then be compared directly to the 
C(Sn, 5B\\) computed from satellite measurements to constrain the nature of the compressible fluctuations. 

Below we present a comparison of four methods to predict and interpret the behavior of C{8n, 5B\\) as 
a function of fa and T: (1) an analytical estimate using MHD eigenfunctions; (2) synthetic spacecraft data 
using MHD eigenfunctions, including contributions from fast, Alfven, and slow waves; (3) synthetic space- 
craft data using Vlasov-Maxwcll eigenfunctions using only the compressible kinetic fast and slow waves; and 
(4) synthetic spacecraft data using Vlasov-Maxwell eigenfunctions, including contributions from fast, Alfven, 
and slow waves. For each of these cases, we explore three combinations of the angular distribution of the 
constituent wave mode power: (a) all isotropically distributed wave power distributions, (b) all critically 
balanced wave power distributions, and (c) isotropic fast wave and critically balanced Alfven and slow wave 
power distributions . Based on intuition gained from numerical simulations of compressible MHD turbu- 



lence (|Cho fc Lazarianl 120031 ). we expect that the third case of isotropic fast waves and critically balanced 
Alfven and slow waves is the most realistic choice; the alternative selections are included to demonstrate the 
sensitivity of the method to different power distributions in wavevector space. 
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All of the synthetic spacecraft data results of C(6n, 6B\\) in this section are relatively insensitive to the 
angle of the probe trajectory with respect to the equilibrium magnetic field, as long as the trajectory is not 
exactly parallel or perpendicular to the field; all results here used an angle of 45° with respect to both the 
equilibrium magnetic field and the x-axis. As a test of the validity of the Taylor hypothesis, we created sets of 
syntehtic data from a time evolving plasma in which the Taylor hypothesis is not assumed. The time evolution 
of each Fourier component for each wave mode was prescribed by its linear frequency and the rate of time 
evolution is parametrized by the ratio of solar wind to Alfven velocities v/va- The correlation C(5n,8Bn) 
from the time evolved and stationary cases were indistinguishable as long as the motion i s sufficiently super- 
Alfve nic (v/va > 10), which is typically satisfied in the solar wind ( Tu fc Marscb 1995 ; Bruno fc Carbonel 



20051 ). Therefore, we assume the Taylor hypothesis is valid and present results derived from stationary 
synthetic data sets. For each choice of /3j and J 7 , the values of C(6n,6B\\) plotted in Figure 0] are the mean 
of 256 ensembles; the error bars are the standard deviation from this statistical averaging procedure. The 
ensembles are used to average over the random phases of the fluctuations. 



4.1. Analytical Estimate of C(Sn,5B») 

To build our intuition about the behavior of C(5n,SB\\) for turbulence modeled as a spectrum of ran- 
domly phased linear wave modes, we use the normalized compressible MHD eigenfunctions presented in 
Appendix lAl to construct an analytical estimator for C(<5n, SB^ ). For a plasma with an equilibrium 

magnetic field Bo = Boz, the density and parallel magnetic field fluctuation for a Fourier mode with wavevec- 
tor k = k sin #x + k cos 9z are generally given by the normalized Fourier coefficients Sn((3, 9) and SBh ((3, 9). 
We can estimate the correlation by integrating the correlations of the Fourier coefficients over the specified 
angular distributions of power for each wave mode, 





6B llm (j3,6) 


d9 


>JiTEm [Sfhn(J3,6)]M) 2 \ 


/(/£ TO 


SB llm (j3,6) 


d9) 2 



where the sum over m includes the contribution from each of the constituent wave modes. The expected 
dependence on (3 is explicit in the expression, while the choice of mode fraction T is implicit in the choice 
of the fraction of power in each of the contributing wave modes. Since the Alfven wave mode in the MHD 
limit has Fourier coefficients Sn = SBu = 0, it is unnecessary to include a contribution from Alfven waves 
to the estimated correlation S(/3). In the top row of Figure HI S(/3) is plotted for the three angular power 
distribution cases: all isotropic (left), all critically balanced (center), and isotropic fast and critically balanced 
slow (right). For each case, we compute the estimator S(/3) for three fractions of fast-to-total compressible 
wave power T = (orange), T = 0.5 (blue), and J- = 1 (black). 

Two qualitative features are immediately apparent from the top row of Figure HI First, the behavior 
for the mixture of fast and slow modes is dominated by the slow mode in the /3 < 1 region and by the fast 
mode in the /3 > 1 region. This result can be understood by comparing the magnitude of the fast and slow 
mode density fluctuations in these regions. In the small j3 region, the slow mode density fluctuations are 
one to two orders of magnitude larger than the fast mode fluctuations; the opposite case holds for the high 
j3 region, where the density fluctuations for the fast mode are much larger than for the slow mode. Second, 
for the case of both isotropic wave power distributions (left), we see a dependence on (3 both for pure fast 
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Fig. 4. — Theoretical predictions of C(Sn, 5B\\ ) using synthetic spacecraft data as a function of ion plasma 
beta fa and fast-to-total compressible wave power T . By row are analytical predictions based on MHD 
eigcnfunctions (first), synthetic spacecraft data using MHD eigcnfunctions (second), synthetic spacecraft 
data using Vlasov-Maxwell eigenfunctions with only fast and slow waves (third), and synthetic spacecraft 
data using Vlasov-Maxwell eigenfunctions with 90% Alfven wave power, and the remaining wave power split 
between fast and slow waves (fourth). By column are presented different angular distributions of wave power: 
all isotropic (first), all critically balanced (second), isotropic fast and critically balanced Alfven and slow 
(third). Results are presented for varying fractions of fast-to-total compressible wave power T: 1.0 (black), 
0.9 (red), 0.75 (green), 0.5 (blue), 0.25 (magenta), 0.1 (grey), 0.0 (orange). The error bars are the standard 
deviation from each ensemble of runs. 
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modes at /3 > 1 and for pure slow modes at /3 < 1. Noting that these pure modes have practically no (3 
dependence for the critical balance case (center), we can surmise that the parallel modes are the cause of the 
deviations seen in the isotropic case. The case of the combination of isotropic fast and critically balanced 
slow modes (right) appears to confirm this finding because the < 1 (slow-wave dominated) region looks 
like the purely critically balanced case, while the /? > 1 (fast-wave dominated) region appears very much like 
the purely isotropic case. These qualitative characteristics of the analytical estimator S(/3) for C(6n,5B^) 
provide a foundation upon which to interpret the synthetic spacecraft data predictions. 



4.2. Synthetic Spacecraft Data Prediction of C(Sn, 5B^) using MHD Fast, Alfven, and Slow 

Eigenmodes 

The compressible MHD cigenfunctions presented in Appendix [5] are used in the procedure outlined in Sj3] 
to generate synthetic spacecraft data to theoretically predict the behavior of C(5n, 8B\\) as a function of the 
synthetic turbulent plasma properties. The synthetic plasma volume is sampled at N = 32 uniformly spaced 
points along a trajectory of length L — vT = ^ x 10 3 pi. Although the MHD Alfven wave has Sh ~ 5B\\ = 0, 
the Alfven wave contribution to the synthetic turbulent plasma is included for completeness, with 90% of the 
turbulent magnetic power given by incompressible Alfven waves. The remaining 10% of the magnetic power 
is split between the MHD fast and slow waves according to the specified fraction of fast-to-total compressible 
wave power T . Tests have shown that the results for C(5n, 6B\\) using MHD cigenfunctions are unaffected by 
the presence or absence of the Alfven wave contribution, as expected. The resulting zero-lag cross correlation 
of the density and parallel magnetic field C(8n,SB\\) from the synthetic MHD plasma is presented in the 
second row of Figure |4] 

We see that C(Sn, 8B\\ ) from the MHD synthetic data displays the same qualitative behavior as the 
analytical estimate S(/3). The dominance in the low (high) (3 regions by the slow (fast) mode behavior is 
evident. A plasma that has 90% of its compressive energy in the fast mode has a slightly negative correlation 
at (3 = 0.1, while a 90% slow mode plasma has a small but positive correlation for /3 = 10. Both of these 
values are drastically different from C(8n,8B\\) at (3 = 1.0. This dominant behavior holds true for all three 
choices of power distributions. In comparing the results from these distributions, we see the marked, and 
expected, lack of (3 dependence for the pure fast and slow modes in the critical balance cases. 



4.3. Synthetic Spacecraft Data Prediction of C(Sn, 5B\\) using only Kinetic Fast and Slow 

Vlasov-Maxwell Eigenmodes 



Unlike the MHD case, the Alfven wave in the Vlasov-Maxwell kinetic theory has a small but non-zero 
fluctuating density and parallel magnetic field component in the incrtial range limit, kpi <C 1. To illuminate 
the contribution of each of the linear kinetic wave modes to C{5n, SB^), in this section we generate a 
synthetic plasma consisting of a spectrum of only kinetic fast and slow wave fluctuations; in the next section, 
the (mostly incompressible) Alfvenic contribution will be included for completeness. 

The sampling of the synthetic plasma volume is the same as for the MHD case, using N — 32 and 



L = ?y x 10 3 pi. Numerical computation of the linear Vlasov-Maxwell eigenfunctions uses the parameters 
" 4 , and Wj/ we 



Tj/Te = 1, = 10 , and Wj/m e = 1836, with B essel function sums evaluated to 100 terms to ensure 
accurate results ( Quataert 1998 ; Howes et al. 20061 ). We have assumed isotropic Maxwellian distribution 
functions for protons and electrons for all Vlasov-Maxwell synthetic spacecraft data results presented in this 
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paper; further exploration of anisotropic temperature distributions may prove fruitful but is beyond the 
scope of the present work. 

The results for C(8n, 8B\\ ) from the Vlasov-Maxwell synthetic plasma with only kinetic fast and slow 
waves are presented in the third row of Figure |4j Comparing to the MHD results (second row), we note an 
apparent similarity in the qualitative behavior but also some noticeable quantitative differences. First, the 
slow mode dominance in the ft < 1 region is still observed for the kinetic plasma, but the ft > 1 region is not 
dominated by the fast mode as in the MHD case. In fact, the mixed modes have very little dependence on 
ft for ft > 1.3. Second, for pure fast and slow modes, the correlation for the isotropic case (left) does have 
a ft dependence in the high and low ft regions respectively, as is seen in the MHD case. However, for the 
pure slow mode, C(5n, 8B\\ ) is never perfectly anti-correlated, reaching a minimum value of ~ —0.9. Third, 
as with the MHD plasma, the pure fast modes are perfectly correlated in the low ft region and become 
drastically less correlated for ft values above the mode conversion (see Appendix [B]) at ft ~ 1.3. 



4.4. Synthetic Spacecraft Data Prediction of C(5n,SB\\) using Kinetic Fast, Alfven, and Slow 

Vlasov-Maxwell Eigenmodes 

In this section, we incorporate an Alfven wave contribution composing 90% of the turbulent magnetic 
power into the turbulent synthetic plasma, with the remaining 10% split between the kinetic fast and slow 
modes. Note that the total compressible wave power used to calculate J- includes only the kinetic fast and 
slow wave contributions and does not include the small contribution to the compressible energy from the 
Alfven wave. Since the density and parallel magnetic field fluctuations of the Alfven wave both have small 
amplitudes, it is expected that the addition of the Alfvenic component will not yield significant quantitative 
changes in C(5n, SB\\). 

Sampling of the synthetic plasma volume and computation of the linear Vlasov-Maxwell eigenfunctions 
is the same as in the previous section. The resulting C(8n, SB\\) for the Vlasov-Maxwell synthetic plasma 
with kinetic fast, Alfven, and slow waves is presented in the fourth row of Figure 21 In comparison to the 
third row of Figure |4j in which the 90% of Alfvenic fluctuation energy is not included, it is clear that the 
inclusion of the dominant Alfvenic component of the turbulence does not lead to significant quantitative 
changes in C(Sn, 5B\\ ). 



4.5. Sensitivity of C(Sn, SB») to Violation of the Taylor Hypothesis 



To test the effect of violating the Taylor hypothesis (|Tavlorl 119381 ) on our results for C(Sn, <5£?||), we 
construct a new time series from a temporally evolving plasma where the Taylor hypothesis is not assumed. 
By abandoning the Taylor hypothesis, we are forced to use equation ([4]) to compute the time series for each of 
the fields computed using our synthetic data method. In this case, each contributing wave mode varies with 
the appropriate linear frequency w m (k) for the particular wavevector k and wave type m. For a particular 
wave mode with plasma-frame frequency w and wavevector k sampled by a moving probe, the normalized, 
Dopplcr-shiftcd wave frequency, given by ZJ' = uj+ (v/va) cos8k v , depends on three dimensionlcss quantities: 
the normalized linear wave frequency uJ = uj/(kvA), the ratio of the probe velocity to the Alfven velocity 
v/va, and the angle Okv between the wavevector of a particular mode k and the probe velocity v, such that 
k • v = kv cos 6kv 
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Fig. 5. — Theoretical prediction of C'(dn, 5B\\) using synthetic spacecraft data as a function of ion plasma 
beta Pi and fast-to-total compressible wave power J- when the Taylor hypothesis is not assumed. Synthetic 
spacecraft data is produced using Vlasov-Maxwcll cigenfunctions with 90% Alfvcn wave power, and the 
remaining wave power split between fast and slow waves. Variation of the ratio of the probe to Alfvcn 
velocity tests the sensitivity of the results to the Taylor hypothesis, where we have tested values v/va = 
oo, 100, 30, 10, 3, 1, where v/va = oo corresponds to the Taylor hypothesis. 



For a synthetic model given by particular spectrum of linear wave modes, each with a given distribution 
of power in wavevector space, once the direction of the probe trajectory is specified, the only remaining 
variable is the ratio of the probe velocity to the Alfven velocity v/va- Therefore, for the same synthetic 
plasma model as specified in £14.41 we may simply vary the value of v/va to observe the effect that violating 
the Taylor hypothesis has on C(5n, SB^ ). The Taylor hypothesis co rresponds to v/va = oo, whereas solar 
wind flows in the near-Earth environment typically have v /v A ^ 10 |Tu fc Marschlll995t bruno fc Carbonel 
20051 ). In Figure [SJ we plot C(8n,6Bn) for values v/va = oo, 100,30, 10,3, 1 for a probe velocity travelling 



at 45° with respect to the mean field direction. For values of v/va > 10, the quantitative effect on the 
correlation is negligible for most of the (/?$, !F) parameter space, with the exception of slight quantitative 
changes for the fast wave dominated cases (J- > 0.9) at > 1. Even for values as low as v/va = E 
the qualitative appearance of the correlation C(6n,5Bu) vs. the ion plasma beta Pi and the fast-to-total 
compressible wave power T is essentially unchanged. Therefore, we conclude that the violation of the Taylor 
hypothesis has little effect on the results of this investigation. 



5. Comparison to Observational Results and Discussion 



Observational constraints suggest that turbulent magnetic power in inertial range turbulence in the solar 
wind consists of approximately 90% Alfven wave s, and the remaining 10% of the power in s ome mixture 
of the compressible kinetic fast and slow waves (iTu fc Marscdll995t Bruno fc Carbonel l2005h . Numerical 
simulations of compressible MHD turbulence suggest that the distribution of turbulent po wer in wavevector 



space is isotropic for the fast waves and critically balanced for the Alfven and slow waves (|Cho fc Lazarian 
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Fig. 6. — Comparison of measured values of the C(Sn, SB\\) cross-correlation (black dots with FWHM error 
bars) to the synthetic data predictions for MHD (left) and kinetic (right) theory given the ratio of fast 
wave to total compressible energy J 7 . Best agreement is with the kinetic theory for J- = 0.00, indicating 
that the compressible component of solar wind turbulence is almost entirely in the kinetic slow mode. The 
compressive behavior is not well described by MHD for any J- '. 



2003). Therefore, we believe the most realistic model of solar wind turbulence using synthetic spacecraft 



data is given by the lower right-hand plot in Figure |4j 

The analysis of the density-parallel magnetic field correl ation using 10 years of Wind spacecraft data is 



discussed in detail in a companion work (jHowes et a.l.ll2012f ) . so we give here only a few brief details. The 



density-parallel magnetic field correlatio n C(5n, SBu) is computed using measurements from the Magn etic 



Field Investigation (MFI) dLeppinJl995l) and the Three Dimensional Plasma (3DP) experiment fLi n 1995) 



on 



the Wind spacecraft in the unperturbed solar wind at 1 AU during the years 1994-2004. Using 300-s intervals 
of ambient solar wind data (corresponding to inertial range scales of approximately kpi € (5 x 10 -3 , 5 x 10~ 2 )), 
the proton density and magnetic field measurements at 3 s cadence are decimated by a factor of 10 (to 30 s 
cadence). Magnetic field measurements are rotated to a field- aligned coordinate system, defined by the local 
mean field direction computed using 100-s windows, to compute SBu , and proton density data 5n = n — tiq 
is detrended over the same time intervals. The zero-lag cross-correlation C(Sn, SB\\) is computed for 119,512 
data intervals. A joint histogram of C(Sn, SB\\) normalized in each /?.; bin is generated, and the peak 
histogram values and FWHM error bars are plotted on top of the theoretical synthetic spacecraft data plots 
of C(5n, 5B\\) from MHD (left) and kinetic (right) theory in Figure |U 

For this direct comparison to the data, the theoretical results from the synthetic spacecraft data are 
computed in the same manner as the spacecraft measurements. In particular, the synthetic plasma volume 
is sampled at only A = 10 uniformly spaced points along a trajectory of length L = vT = x 10 2 pi, which 
corresponds to a range of scales kpi £ (4.8 x 10~ 3 , 4.8 x 10~ 2 ). The smaller number of timesteps and shorter 
total sampled length leads to a larger standard deviation using 256 ensembles than the Vlasov-Maxwell 
results in fj4]and to slight quantitative changes in the C(8n,6B\\) curves. 
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The agreement between the C(8n, 8B\\ ) from Wind spacecraft data and synthetic data curve for F = 
in Figure [5] is striking, indicating that the observed correlation is consiste nt with a statistica lly negligible 
kinetic fast wave energy contribution for the large sample used in this study (jHowes et al.ll2012l) . We do note 
that a very small fraction of the intervals have C (8n, 6B \\ ) > 0, poss ibly indicating a minority population of 



fast waves. As discussed in our companion paper (jHowes et al.ll2012l ). this result has important consequences 
for the turbulent cascade of energy from large to small scales: since only the fast wave turbulent cascade is 
expected to nonlinearly transfer energy to whistler waves at kpi > 1, and the frequency mismatch between 
the fast and either the Alfven or slow waves should prevent non-linear coupling, our analysis suggests that 
there is little or no transfer of large scale turbulent energy through the incrtial range down to whistler waves 
at small scales. 

Comparison of the observationally measured C{8n,8B\\) in Figure |6] with both the MHD and kinetic 
results leads to another important conclusion of this study: the nature of the compressible fluctuations in 
the solar wind inertial range is not well modeled by MHD theory. This result is not surprising given that 
the solar wind is weakly collisional at inertial range scales, so a kinetic description is necessary to model 
accurately the compressible fluctuations. 

Our study also enables us to constrain the wavevector distribution of slow wave energy by comparing the 
observational results with the three plots in the bottom row of Figure UJ The leftmost plot has an isotropic 
distribution of slow wave power and shows two features in the T = curve not seen in the observed data: (1) 
a bump at ~ 1 and (2) a slight increase in the correlation value at < 0.3. These quantitative changes 
in the T = curve appear the same with or without the Alfven waves at Pi < 1 (compare the left plots in 
the third and fourth row of Figure |4j), so this effect is not due to the isotropic distribution of Alfven waves. 
Therefore, the slow waves appear to be anisotropically distri buted, and may be well described by critical 
balan ce as suggested by theories of slow wave passive advection ( Maron fc Goldreichl2001 ; Schekochihin et al 



20091 ) . However, the uncertainty of the Wind measurements cannot definitively rule out the possibility of an 
isotropic distribution for the slow mode. 



One of the key points of our companion work (jHowes et al.l l2012j ) is the interpretation that the com- 
pressible fluctuations in the solar wind consist of a critically balanced distribution of kinetic slow wave 
fluctuations. Previous analyses have generally dismissed the possibility of kinetic slow waves because, in an 
isotropic Maxw ellian plasma with warm ions, the collisionless damping via free-streaming along the magnetic 
field is strong ( Barnes! 1966 ). However, the damping rate of the slow waves for the nearly perpendicular 
wavevectors of a critically balanced distribution is proportional to the parallel component of the wavevector, 
7 oc fen. This feature can be seen in the right column of Figure [T] where the slow wave damping rate (blue) 
is zero at perpendicular wavevectors 9 = 90° (open circles) and increases as the wavevector angle decreases. 
For exactly perpendicular wavevectors, the damping rate drops to zero — this perpendicular limit of the slow 
wave correspon ds to an undamped, non-p ropagating pressure-balanced stru cture (PBS). It has been de rived 
theoretically by I Schekochihin et al.l (120091 ) and demonstrated numerically bv lMaron fc Goldreichl (j200ll) that 
the Alfven wave dynamics advects and cascades the slow waves, so the energy cascade rate of the slow waves 
is related not to the slow wave frequency, but to the Alfven wave frequency. Therefore, although the slow 
wave fluctuations at the high fc|| boundary of critical balanced distribution (along fen ~ k^^k 2 ^ 3 ) may suffer 
strong collisionless damping, the more nearly perpendicular slow waves may be cascaded to smaller scales 
on the timescale of the Alfvenic turbulence, while the collisionless damping of these modes remains weak. 
This could lead to a slow wave energy distribution that is more anisotropic than critical balance, possibly 
with fen <C feg k i 3 . 



Using synthetic spacecraft data, we can test whether this picture of a distribution of slow wave power 
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Fig. 7. — Comparison of synthetic data predictions for C(Sn, 6B\\) with isotropic fast mode power, critically 
balanced Alfven mode power, and either critically balanced (solid) or more anisotropic than critical balance 
(dashed) slow mode power. For the more anisotropic than critical balance distribution, all modes with 

2 1/3 

kit > k fc_L are set to zero, where fco is the minimum wavenumbcr of the wavenumber domain. The 
difference between the lines is attributable to our choice for specifying the fraction of wave power individually 
for each Fourier component: the more anisotropic distribution contains less total slow wave power. The lines 
for the J- = are nearly identical for both power distributions, signifying that C(5n,5Bu) can not be used 
to differentiate between critical balance and more anisotropic distributions. 



that is more anisotropic than critical balance is consistent with the observational findings in Figure |6] To 
do this, we repeat the analysis in the lower right plot of Figure |4] but with a more anisotropic distribution of 
slow waves given by ku < fc^ 3 ^ 3 . The results of this test, shown in Figure [71 demonstrate that the J- = 
curve does not change significantly for a more anisotropic slow wave distribution. Therefore, the picture of 
the solar wind compressible fluctuations comprised of slow wave fluctuations that are more anisotropically 
distributed than critical balance is consistent with the observational data. 

One of the limitations of this investigation is that the Vlasov-Maxwell eigenfunctions are computed 
assuming isotropic Maxellian equilibrium distributions. Plasma measurements in t he solar wind freque ntly 



find anisotropic temperature distributions and non-Maxwellian tails at high energy (jMarsch 



1991 



20061 ). In 



particular, it is possible that mirror modes may contribute to the measured anti-correlation between the 
density and parallel magnetic field, although an analysis of 10 years of Wind data show that only a small 
fraction o f the solar wind measurements occupy the region of parameter space near the mirror instability 
threshold ( Bale et aDl2009h . How the density-parallel magnetic field correlation is altered by these conditions 
is an important issue to be addressed in future work. 
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6. Conclusion 

In this paper, we have presented the first application of kinetic plasma theory to analyze and interpret 
the compressible fluctuations in the inertial range of solar wind turbulence. This novel approach is motivated 
by the fact that the dynamics in the solar wind plasma is weakly collisional, a limit in which MHD theory, 
used in all previous related investigations, is formally invalid. Investigation of the compressible fast and slow 
wave modes requires a kinetic description to resolve both the wave dynamics and the collisionless kinetic 
damping mechanisms. 

We identify quantitatively the linear kinetic wave modes of Vlasov-Maxwell theory that correspond to 
the linear MHD fast and slow wave modes, and verify that the general qualitative properties of the density- 
parallel magnetic field correlation C(Sn, 6B\\) — that fast waves are positively correlated, and slow waves are 
negatively correlated — remain unchanged for the weakly collisional conditions of the solar wind plasma. 

We then describe the procedure used to generate synthetic spacecraft data used to interpret actual 
single-point spacecraft measurements of turbulence in the solar wind. We define and discuss the quasi- 
linear premise — that some properties of magnetized plasma turbulence can be understood by modeling the 
turbulence as a collection of randomly phased, linear waves — upon which the method of synthetic spacecraft 
data is based. Theoretical arguments for the validity of the quasi-linear premise are presented ; a review of 



suppo rting and conflicting evidence for the quasi- linear premise in the literature is presented in lHowes et al 



(|201ll ). We outline how the synthetic plasma data cubes are used to generate synthetic reduced time series 
of single-point measurements that can be analyzed using the same procedures as the actual spacecraft 
measurements. Comparison of the results of the analyses of both the synthetic a nd actual spacecra ft data 



enables novel physical interpretations of the solar wind turbulence measurements (IHowes et al.ll2012l) 



Next, we use the synthetic spacecraft data method to predict the characteristics of the density-parallel 
magnetic field correlation C(dn, SBu) as a function of plasma /3j in the presence of a turbulent spectrum of 
wave modes. Comparison of these synthetic predictions of C(6n, SB^), based on both linear MHD and linear 
Vlasov-Maxwell eigenf unctions, to the results from an analysis of a 10-year data set of observations from the 
Wind spacecraft leads to the following conclusions: 

1. The predicted C(<5n, 5B\\) and its dependence on plasma fii using linear MHD eigenfunctions is signif- 
icantly different from the prediction using linear Vlasov-Maxwell eigenfunctions. Only the prediction 
based on kinetic theory appears to agree with the spacecraft measurements, leading to the expected 
conclusion that MHD theory is inadequate to describe the compressible fluctuations in the weakly 
collisional solar wind. 

2. Strong a posteriori evidence for the validity of the quasi-linear premise is provided by the striking 
agreement between the observationally determined C(Sn, 5B\\) over a very large statistical sample and 
the predicted C(5n, SB\\) based on synthetic spacecraft data. 



The observed C(8n,5B\\) computed in a companion work (jHowes et al.l 120121 ) is consistent with a 
statistically negligible kinetic fast wave energy contribution for the large sample used in this study. 
Note, however, that the our companion work also found that a very s mall fraction of the intervals have 



C(5n, 5B\\) > 0, possibly indicating a trace population of fast waves (jHowes et al.ll2012l ). 



4. The quantitative dependence of C{5n, 5B\\ ) on the ion plasma beta fit provides evidence that the 
slow wave fluctuations are not isotropically distributed, but rather have an anisotropic distribution, 
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that is possibly given by the condition of critical balance (IGoldreich fc Sridharlll995l ) or that is more 
anisotropic than critical balance. 



In conclusion, our analysis using kinetic theory to interpret the compressible fluctuations motivates the 
following physical model of the turbulent fluctuations in the solar wind incrtial range. In this model, the 
solar wind inertial range consists of a mixture of turbulent fluctuations, with 90% of the energy due to 
incompressible Alfvenic fluctuations and the remaining 10% of the energy due to compressible slow wave 
fluctuations. The Alfvenic turbulent power is distributed anisotropically in wave vector space according to 
critical balance. The turbulent Alfven wave dynamics advects and cascades the slow wave fluctuations to 
smaller scales at the Alfven wave frequency. The slow wave turbulent power may be either critically balanced, 
or more anisotropic than critical balance due to collisionless damping of the slow wave fluctuations. 

Since only the fast wave turbulent cascade is expected to nonlinearly transfer energy to whistler waves 
at kpi > 1, and the frequency mismatch between the fast and either the Alfven or slow waves should prevent 
non-linear coupling, there is little or no transfer of large scale turbulent energy through the inertial range 
down to whistler waves at small scales. Therefore, any whistler wave fluctuations at s cales kpj > 1 must 



be generated by some other process, e.g., kinetic temperature anisotropy instabilities (jKasper et al.l 12002 ; 



Hellinger et"alll2006t 



Bale et al 



2009 1) or k i netic drift instabilities driven by differential flow between protons 



and alpha particles (jMcKenzie et al.lll993t iKasper et al.ll2008t iBourouaine et al]|2011l ). 



Finally, the lack of statistically significant fast wave energy has important implications for efficient nu- 
merical modeling of solar wind turbulent fluctuations. This work demonstrates clearly the importance of a 
kinetic approach to model adequately the turbulent fluctuations, yet a general kinetic numerical treatment — 
e.g., the Particle-In-Ccll (PIC) method — in three spatial dimensions (required for phys ically relevant model- 
ing of the dominant nonlinear interactions in solar wind turbulence (jHowes et al.ll201ll )) is too computation- 
ally costly to be presently feasible. Fortunately, it is possible to perform kinetic numerical simulations of solar 
wind turbulence in three spatial dimensions using gyrokinetics, a rigorous, low-frequency, anisotropic limit 
of kin etic theory ( Rutherford &: Frieman 1968 : Frieman &; Chen 1982 : Howes et al. 20061 : Schekochihin et al 



20091). In the derivation of the gyrokinetic equation, the crucial step is an averaging over the particle gy- 
rophase, which leads to a theory with the following properties: the fast/whistler wave and the cyclotron 
resonances are discarded; all finite Larmor radius effects and collisionless dissipation via the Landau reso- 
nance are retained; and one of the dimensions of velocity in phase space is eliminated, reducing the particle 
distribution function from six to five dimensions. It has been previously pointed out that one cannot rule 
out the contribution of fast wave or whistler wave physics t o solar wind turbulenc e, and that that therefore 
gyrokinetics is an incomplete description of the turbulence ( Matthaeus et al. 20081) . The novel observational 
analysis presented here suggests that the fast wave, in fact, does not play a statistically significant role in 
the turbulent dynamics of the incrtial range, and that therefore a gyrokinetic approach sufficiently describes 
all important physical mechanisms in the solar wind inertial range. 



K. G. K. and G. G. H. thank Jack Scudder for inspiring discussions. The work has been supported by 
NSF CAREER Award AGS-1054061 and NASA NNX10AC91G. 
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A. Normalized MHD Linear Eigenfunctions 



The ideal, compressible MHD equations are given by the continuity equation, 

^ + U • Vp = -pV ■ u, 

the momentum equation, 

(B • V)B 



D- 



4tt 



the induction equation, 



and an adiabatic equation of state, 



- = Vx(UxB), 



dv 

-jf + U • Vp = -7pV • U. 
ot 



(Al) 



(A2) 



(A3) 



(A4) 



Here the MHD fluid is completely described by its mass density p = n(m,i + m e ), fluid velocity U, magnetic 
field B, and scalar thermal pressure p. The adiabatic index for a proton and electron plasma is 7 = 5/3. 
This system of equations can be derived rigorously from plasma kinetic theory in the MHD limit of strong 
collisionality v 3> w, large scales compared to the ion Larmor radius kpi <C 1, and non-relativistic conditions 
dKulsrudlll98i . 



Although the weakly collisional conditions of the solar wind plasma violate the strong collisionality 
formally required for the validity of MHD, this system has nevertheless been widely used to study the 
turbulent dynamics of the solar wind inertial range. That this simplified approach has not met with 
widespread failure is likely due to the fact that the turbulent dynamics in the inertial range is dominated by 
Alfvenic motions (jBelcher fc Davislll97ll : iTu fc Marschlll995t iBruno fc Carbonell2005h . Alfvenic fluctuations 
at inertial range scales are incompressible, and it has been shown that for anisot ropic fluctuations with 
<S! fcj_ (the anisotropy generally observed in the magnetized plasma turbulence ( Robinson &: Rusbridgel 
197ll ; IZweben et alll979t iMontgomerv fc TurneJll98ll: ISahraoui et~all201fi Ichen et alJboill )). the dvnam- 
ics of Alfvenic fluctua tions are rigorously descr ibed by the equations of reduced MHD, even under weakly 
collisional conditions (jSchekochihin et al.1 120091 ) . The compressible fast and slow wave fluctuations, on the 
other hand, are subs tantially modified in weakly collisional conditions, so their dynamics must be treated 
using kinetic theory ( Kulsrud 1983 ; Schekochihin et al. 20091) . as has been demonstrated in the investigation 
presented here. 

In this appendix, we derive the compressible MHD linear dispersion relation and eigenfunctions in 
dimensionless units constructed specifically for our study of compressible solar wind fluctutations. Without 
loss of generality, we specify a wavevector k = fc^x + fc||Z and an equilibrium magnetic field Bo = Bqz. We 
separate mean from fluctuating quantities using p — po + Sp, U = 5XJ, B = Bo + <5B, and p = pg + 5p, 
where we specify that there is no mean fluid velocity. We then linearize equations (|A1|) (| A4|) and Fourier 
transform the equations in both space and time. Next, we convert each of the variables to dimensionless 
units and define the dimensionless frequency and other parameters as follows: 



5n = 



Sp 

sXj 

5B 

5p 



Sp/po 
SV/va 
dB/B 



ft 



Uj/{kVA) 

IPolPo 
Bl/inpo 



(A5) 



-24- 



In this normalization, the linear dispersion relation takes the form 



- 2 — 2 9) [w 4 - uJ 2 (l + P) + p cos 2 9]=0 (A6) 



[UJ — cos 



where the first factor in parentheses corresponds to the two Alfven wave solutions, and the second factor in 
the brackets corresponds to the two slow and two fast wave solution^. This demonstrates the important 
property that the normalized compressible MHD linear dispersion relation depends on only two parameters, 
Lo/ikvA) = wmhd(/3,0): the plasma beta /3 and the angle 9 between the local mean magnetic field Bo and 
the direction of the wavevector k. 

The linear Alfven wave solutions have frequency 

uj 2 = cos 2 9 (A7) 
and eigenfunctions specified in terms of SU y given by 

JWy = ±SU^ (A8) 
7m = WT X = W\\ = W x = 5B^ = fy = (A9) 

The linear slow and fast wave solutions have frequency 



_ 2 1 + ± y/(l + /?) 2 - 4/W 

oj = , (A10) 



Sn = -» ^— (All) 



where the plus sign corresponds to the fast wave, and the minus sign to the slow wave. The eigenfunctions, 
specified in terms of 8Bu , are given by 

w 2 - p cos 2 6 

JB~ X = -W^cot9 (A12) 

=lJ5B^ csc 9 (A13) 

— — uPcos9 — — / A1 ^ 

5U\\ = — 8B\\ ( Al4 ) 

" uj 2 -Pcos 2 9 " V ; 

5Uy = SB^ = (A15) 

Since the synthetic data sets are created by specifying a spectrum of fluctuations over some range 
of angles, the statistical correlations of the fluctuations depend on only three factors: (a) angular power 
distribution, (b) wave mode fraction, and (c) plasma /3. 



B. Determining Kinetic Counterparts of Slow vs. Fast MHD Waves 



The kinetic fast and slow modes are both part of the double- valued magnctosonic solution (jKrauss-Varban et al 



1994T ) and are separated by a branch cut from 9 w 30° to 0° at /3j w 1.3. This branch cut is similar to the 



magnetic and acoustic transition for the MHD fast and slow modes at j3 = 1 for near parallel wavevectors. 



5 Note that entropy mode solution uj = has already been removed from this dispersion relation since our focus here is on 
the propagating linear wave modes. 
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The kinetic fast and slow modes behave in an analogous fashion: the near parallel kinetic fast mode has the 
same phase velocity as the Alfven mode for ft < 1.3; and the near parallel kinetic slow mode matches the 
Alfven mode for ft > 1.3 (see the center column of Figure [1]). In this work, we label wave modes using a 
scheme that identifies the modes at 6 = 90°, where there is no ambiguity in identification, and then follows 
the linear dispersion relation as 6 is decreased to zero. 

The compressible inertial range linear wave modes with near parallel wavevectors can be characterized 
as being either magnetic or acoustic. The magnetic modes have phase velocities similar to the Alfven wave, 
while the acoustic modes are more heavily damped than either the Alfven or the magnetic modes. From this 
point of view, the fast mode is the magnetic mode and the slow mode is the acoustic mode for ft < 1.3, while 
the converse is true for ft > 1.3. Interestingly, for ft < 1.3, the magnetic (fast) mode has C(5n,8B\\) ~ 1 
and the acoustic (slow) mode has C(5n, 6B\\) ~ —1; but, for ft > 1.3, both modes have C(5n, 5B\\) ~ (see 
Figure [2]) . It is this transition that is responsible for some of the qualitative features in the lower two rows of 
Figure SI In particular, the abrupt jump on the J- = 1 curve (black) at ft ~ 1.3 in the left and right columns 
(where fast waves are isotropically distributed and so include this region of wavevector space) and in the 
J- = curve (orange) at ft ~ 1.3 in the left column (where the slow waves are isotropically distributed) are 
due to this transition. 

This identification of modes has important physical ramifications, especially at scales near the transition 
to the dissipation range. As length scales decrease towards the ion inertial length, it is the magnetic mode, 
and not strictly the kinetic fast mode, that transitions to the parallel whistler wave. However, the results 
presented in this paper suggest a lack of parallel slow wave energy as well as statistically negligible fast wave 
energy in the solar wind inertial range, so the difficulties of mode identification in this region of parameter 
space may not be problematic for studies of solar wind turbulence. 

REFERENCES 

Bale, S. D., Kasper, J. C, Howes, G. G., Quataert, E., Salem, C, & Sundkvist, D. 2009, Phys. Rev. Lett., 
103, 211101 

Barnes, A. 1966, Phys. Fluids, 9, 1483 

Baumjohann, W., & Treumann, R. A. 1996, Basic space plasma physics (London: Imperial College Press) 
Bavassano, B., Pietropaolo, E., & Bruno, R. 2004, Annales Geophysicae, 22, 689 
Belcher, J. W., & Davis, L. 1971, J. Geophys. Res., 76, 3534 
Boldyrev, S. 2006, Phys. Rev. Lett., 96, 115002 

Bourouaine, S., Marsch, E., & Neubauer, F. M. 2011, Astrophys. J. Lett., 728, L3+ 

Bruno, R., & Carbone, V. 2005, Living Reviews in Solar Physics, 2, 4 

Burlaga, L. F. 1968, Sol. Phys., 4, 67 

Burlaga, L. F., & Ogilvie, K. W. 1970, Sol. Phys., 15, 61 

Burlaga, L. F., Scudder, J. D., Klein, L. W., & Isenberg, P. A. 1990, J. Geophys. Res., 95, 2229 



-26 - 



Chen, C. H. K., Mallet, A., Yousef, T. A., Schekochihin, A. A., & Horbury, T. S. 2011, Mon. Not. Roy. As- 
tron. Soc., 415, 3219 

Cho, J., & Lazarian, A. 2003, Mon. Not. Roy. Astron. Soc., 345, 325 

Cho, J., & Vishniac, E. T. 2000, Astrophys. J., 539, 273 

Coleman, Jr., P. J. 1968, Astrophys. J., 153, 371 

Dmitruk, P., & Matthaeus, W. H. 2009, Phys. Plasmas, 16, 062304 

Elsasser, W. M. 1950, Physical Review, 79, 183 

Forman, M. A., Wicks, R. T., & Horbury, T. S. 2011, Astrophys. J., 733, 76 
Frieman, E. A., & Chen, L. 1982, Phys. Fluids, 25, 502 

Galtier, S., Nazarenko, S. V., Newell, A. C, & Pouquet, A. 2000, J. Plasma Phys., 63, 447 
Goldrcich, P., & Sridhar, S. 1995, Astrophys. J., 438, 763 
— . 1997, Astrophys. J., 485, 680 

Hellinger, P., Travnicek, P., Kasper, J. C, & Lazarus, A. J. 2006, Geophys. Res. Lett., 33, 9101 
Higdon, J. C. 1984, Astrophys. J., 285, 109 

Horbury, T. S., Forman, M., & Oughton, S. 2008, Phys. Rev. Lett., 101, 175005 
Howes, G. G. 2009, Nonlin. Proc. Geophys., 16, 219 

Howes, G. G., Bale, S. D., Klein, K. G., Chen, C. H. K.and Salem, C. S., & TenBarge, J. M. 2012, 
Astrophys. J. Lett., submitted 

Howes, G. G., Cowley, S. C, Dorland, W., Hammett, G. W., Quataert, E., & Schekochihin, A. A. 2006, 
Astrophys. J., 651, 590 

— . 2008a, J. Geophys. Res., 113, A05103 

Howes, G. G., Dorland, W., Cowley, S. C, Hammett, G. W., Quataert, E., Schekochihin, A. A., & Tatsuno, 
T. 2008b, Phys. Rev. Lett., 100, 065004 

Howes, G. G., Tenbarge, J. M., & Dorland, W. 2011, Phys. Plasmas, 18, 102305 

Howes, G. G., TenBarge, J. M., Dorland, W., Quataert, E., Schekochihin, A. A., Numata, R., & Tatsuno, 
T. 2011, Phys. Rev. Lett., 107, 035004 

Howes, G. G., TenBarge, J. M., & Klein, K. G. 2011, Astrophys. J., in preparation 

Hunana, P., Laveder, D., Passot, T., Sulem, P. L., & Borgogno, D. 2011, Astrophys. J., 743, 128 

Kasper, J. C, Lazarus, A. J., & Gary, S. P. 2002, Geophys. Res. Lett., 29, 20 

— . 2008, Phys. Rev. Lett., 101, 261103 

Kellogg, P. J., & Horbury, T. S. 2005, Ann. Geophys., 23, 3765 



-27- 



Krauss-Varban, D., Omidi, N., & Quest, K. B. 1994, J. Geophys. Res., 99, 5987 

Kulsrud, R. M. 1983, in Handbook of Plasma Physics, Vol. 1, Basic Plasma Physics I, ed. A. A. Galeev & 
R. N. Sudan (North Holland), 115-145 

Lepping, R. 1995, Space Sci Rev, 71, 207 

Lin, R. P. 1995, Space Sci Rev, 71, 125 

Lithwick, Y., & Goldreich, P. 2003, Astrophys. J., 582, 1220 
Luo, Q. Y., & Wu, D. J. 2010, Astrophys. J. Lett., 714, L138 
Maron, J., & Goldreich, P. 2001, Astrophys. J., 554, 1175 

Marsch, E. 1991, in Physics of the Inner Heliosphere II. Particles, Waves and Turbulence., ed. E. Schwenn, 
R. andMarsch (Berlin: Springer- Vcrlag) , 45-133 

Marsch, E. 2006, Living Rev. Solar Phys., 3 

Matthacus, W. H., Klein, L. W., Ghosh, S., & Brown, M. R. 1991, J. Geophys. Res., 96, 5421 
Matthaeus, W. H., Servidio, S., & Dmitruk, P. 2008, Phys. Rev. Lett., 101, 149501 

McComas, D. J., Barraclough, B. L., Gosling, J. T., Hammond, C. M., Phillips, J. L., Neugebauer, M., 
Balogh, A., & Forsyth, R. J. 1995, J. Geophys. Res., 100, 19893 

McKenzie, J. F., Marsch, E., Baumgaertel, K., & Sauer, K. 1993, Ann. Geophys., 11, 341 

Montgomery, D., Brown, M. R., & Matthaeus, W. H. 1987, J. Geophys. Res., 92, 282 

Montgomery, D., & Matthaeus, W. H. 1995, Astrophys. J., 447, 706 

Montgomery, D., & Turner, L. 1981, Phys. Fluids, 24, 825 

Narita, Y., Gary, S. P., Saito, S., Glassmeier, K.-H., & Motschmann, U. 2011, Geophys. Res. Lett., 38, 
L05101 

Ng, C. S., & Bhattacharjee, A. 1996, Astrophys. J., 465, 845 
— . 1997, Phys. Plasmas, 4, 605 

Parashar, T. N., Servidio, S., Breech, B., Shay, M. A., & Matthaeus, W. H. 2010, Phys. Plasmas, 17, 102304 
Podesta, J. J. 2009, Astrophys. J., 698, 986 
Quataert, E. 1998, Astrophys. J., 500, 978 

Reisenfeld, D. B., McComas, D. J., & Steinberg, J. T. 1999, Geophys. Res. Lett., 26, 1805 
Roberts, D. A. 1990, J. Geophys. Res., 95, 1087 

Roberts, D. A., Goldstein, M. L., Klein, L. W., & Matthaeus, W. H. 1987a, J. Geophys. Res., 92, 12023 
— . 1987b, J. Geophys. Res., 92, 11021 

Robinson, D. C, & Rusbridge, M. G. 1971, Phys. Fluids, 14, 2499 



-28 - 



Rutherford, P. H., & Frieman, E. A. 1968, Phys. Fluids, 11, 569 

Sahraoui, F., Goldstein, M. L., Belmont, G., Canu, P., & Rezcau, L. 2010, Phys. Rev. Lett., 105, 131101 

Schekochihin, A. A., Cowley, S. C, Dorland, W., Hammett, G. W., Howes, G. G., Quataert, E., & Tatsuno, 
T. 2009, Astrophys. J. Supp., 182, 310 

Shebalin, J. V., Matthaeus, W. H., & Montgomery, D. 1983, J. Plasma Phys., 29, 525 

Sridhar, S., & Goldreich, P. 1994, Astrophys. J., 433, 612 

Stix, T. H. 1992, Waves in Plasmas (New York: American Institute of Physics) 

Taylor, G. I. 1938, Proc. Roy. Soc. A, 164, 476 

Tu, C.-Y., & Marsch, E. 1994, J. Gcophys. Res., 99, 21481 

— . 1995, Space Sci. Rev., 73, 1 

Vellante, M., & Lazarus, A. J. 1987, J. Geophys. Res., 92, 9893 

Wicks, R. T., Horbury, T. S., Chen, C. H. K., & Schekochihin, A. A. 2010, Mon. Not. Roy. Astron. Soc, 
407, L31 

Yao, S., He, J.-S., Marsch, E., Tu, C.-Y., Pedersen, A., Reme, H, & Trotignon, J. G. 2011, Astrophys. J., 
728, 146 

Zweben, S. J., Menyuk, C. R., & Taylor, R. J. 1979, Phys. Rev. Lett., 42, 1270 



This preprint was prepared with the AAS IATjTjX macros v5.2. 



